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Abstract. High resolution {R ~ 110 000), very high S/N 
spectra centered on the 6707.8 A Lii line have been ob- 
tained with the ESO Coude Echelle Spectrometer for five, 
metal-poor (—0.8 < [Fe/H] < —0.6) disk stars in the 
turnoff region of the HR-diagram. The instrumental and 
stellar atmospheric line broadening have been determined 
from two unblended iron lines and used in a model atmo- 
sphere synthesis of the profile of the Li i line as a function 
of the lithium isotope ratio. This has led to a detection 
of ^Li in HD 68284 and HD 130551 with ^Li/^Li ~ 0.05, 
whereas the other stars, HR2883, HR3578 and HR8181, 
have ^Li/^Li close to zero. 

By comparing Tcg-values and absolute magnitudes 
based on Hipparcos parallaxes with recent stellar evolu- 
tionary tracks, the masses of the stars have been derived. 
It is shown that the two stars with ^Li present have a 
significantly higher mass, M/AAq ~ 1.05, than the other 
three stars for which values between 0.85 and 1.0 J^q are 
obtained. 

The results are discussed in terms of models for the 
galactic evolution of the light elements and depletion of 
the lithium isotopes in stellar envelopes. It is shown that 
the measured ^Li abundances are in agreement with stan- 
dard cosmic ray production of ^Li in the galactic disk and 
a moderate depletion (0.5 dex) in the stars. Recent models 
for the evolution of ^Li including a + a fusion reactions 
and predicting a high lithium isotopic ratio, ^Li/^Li ~ 0.3 
at [Fe/H] — —0.6, require a high degree of ^Li depletion 
(~ 1.0 dex) to fit the observations. Furthermore, these 
models imply a ^Li abundance about 0.2 dex higher than 
observed for metal-poor disk stars. 

Key words: Stars: abundances - Stars: atmospheres - 
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1. Introduction 

Measurements of the abundance of the ^Li isotope in stel- 
lar atmospheres are of considerable interest and have at- 
tracted much attention since the first detection of ^Li 
in the metal-poor turnoff star HD 84937 by Smith et al. 
( 1993| ). The reason for this interest is threefold: 

i) Detection of ®Li in halo turnoff stars puts strong 
limits on the possible depletion of ^Li, and thus allows bet- 
ter determination of the primordial '^Li abundance from 
the observed Li abundance of stars on the 'Spite plateau'. 



(Copi et al. |1997| , Pinsonneault et al. |1998D 

ii) ^Li abundances as a function of [Fe/H] provide an 
additional test of theories for the production of the light 
elements Li, Be and B by interactions between fast nuclei 



and ambient ones (Ramaty et al. 1996, Yoshii et al. 1997, 
Fields & Olive |l999a"^ Vangioni-Flam et al. |1999[ ). 

Hi) Information on depletion of ^Li as a function of 
stellar mass and metallicity puts new constraints on stellar 
models in addition to those set by ''Li depletion. This is 
so because the proton capture cross section of ^Li is much 
larger than that of ^Li. Hence, at a given metallicity there 
will be a mass interval, where ^Li but not ^Li is being 
destroyed according to standard stellar models (Chaboyer 



1994) 



Altogether, ^Li abundances may contribute to the 
study of such different fields as Big Bang nucleosynthesis, 
cosmic ray physics and stellar structure. It will, however, 
require a rather large data set of ^Li abundances to get 
information in all these areas. The most metal-poor stars 
around the turnoff are of particular interest in connection 
with the determination of the primordial ''Li abundance, 
whereas more metal-rich halo stars and disk stars are of 
interest for the study of the formation and astration of the 
light elements. 

Recent studies of ^Li abundances have concentrated on 
halo stars. Following the first detection of ^Li in HD 84937 



by Smith et al. (1993) at a level corresponding to an iso- 
topic ratio of ^Li/^Li ~ 0.05, Hobbs & Thorburn ( |1994 



1997) have confirmed the detection, and found upper lim- 



i ts of ^Li/''Li for 10 stars. More recently. Smith et al. 
( 1998 ) report the probable detection of ^Li in another 
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halo star BD +26 3578 with about the same metal abun- 
dance, mass and evolutionary stage as HD 84937, and give 
tight upper limits of ^Li/^Li for 7 additional stars. Fi- 



nally, Cayrel et al. (1999a) have observed the Lii line in 
HD 84937 with very high S/N and confirmed the presence 
of ^Li with a high degree of confidence. 

In the case of disk stars there has not been any sys- 
tematic search for ^Li since the studies of Andersen et al. 
( |1984D and Maurice et al. (1984). In these papers an upper 
limit of ^Li/^Li of about 0.10 is set for about 10 disk stars 
ranging in metallicity from —1.0 to -1-0.3. The meteoritic 
^Li/^Li ratio is close to 0.08 (Anders & Grevesse [1989[ ) 
and the interstellar ratio is similar - possibly with signifi- 



cant variations (Lemoine et al. [1995 ). For metal-poor disk 
stars the ratio may be considerably higher than in the 
solar system. According to recent models for the ga lactic 
evolution of th e light ele ments (Vangioni-Flam et al. 1999, 
Fields & Olive 1999a, b ) the ^Li/^Li ratio reaches a max- 
imum of about 0.3 at a metallicity of [Fe/H] ~ —0.5. At 
higher metallicities the ratio decreases due to the produc- 
tion of ''Li in AGB stars, novae and su pernovae of type II 
by th e z/-process (Matteuc ci et a l. 1995 , Woosley & Weaver 
1995 , Vangioni-Flam et al. 1996| ). Hence, it seems well jus- 



tified to look for ^Li in the metal-poor disk stars. Any de- 
tection will provide important constraints of the chemical 
evolutionary models, and with a large set of data it may 
also be possible to constrain the degree of ^Li depletion 
as a function of stellar mass and metallicity. 

In the present paper we present results for the ^Li/^Li 
ratio for five metal-poor disk stars ranging in metallicity 
from -0.8 to -0.6. Very high S/N spectra of the Li i 6708A 
resonance line are presented in Sect. 3 and analyzed with 
model atmosphere techniques in Sect. 4. This has led to 
a rather clear detection of ^Li in the two stars with the 
highest masses and tight upper limits for ^Li in the other 
stars. The consequences of these results are discussed in 
Sect. 5. 

2. Selection of stars 

The program stars were selected from the large survey of 



nearby disk dwarfs by Edvardsson et al. (1993), which in- 
cludes accurate values of atmospheric parameters, abun- 
dances, kinematics and ages of 189 main sequence stars 
distributed in metallicity from [Fe/H] = —1.0 to -1-0.3. 
As we wanted to avoid stars formed from interstellar gas 
greatly enriched in ^Li from AGB stars or novae, the first 
condition for including a star was [Fe/H] ^ —0.5. Next, 
only stars with T^s ^ 5900 K were included in order to 
maximize the chance of survival of some ^Li. Finally, the 
program had to be limited to a few of the brightest stars 
in order to be able to reach the high S/N that is needed 
to determine the lithium isotope ratio from the profile of 
the 6707.8 A Li I line. 

The five stars observed are listed in Table 1. The ef- 
fective temperature (derived from the b — y color index) , 



Table 1. List of program stars including the visual mag- 
nitude mv and atmospheric parameters adopted from Ed- 
vardsson et al. ( 1993D 



HR 



HD 



mv 



[K] 



logp [Fe/H] 



[kms 



rmcro 
ll 



2883 


59984 


5.93 


5976 


4.18 


-0.75 


1.7 


3578 


76932 


5.80 


5965 


4.37 


-0.82 


1.4 


8181 


203608 


4.23 


6139 


4.34 


-0.67 


1.6 




68284 


7.77 


5883 


3.96 


-0.59 


1.9 




130551 


7.16 


6237 


4.25 


-0.62 


1.8 



the logarithmic surface gravity, the iron abundance, and 
the microturbulence velocity, are taken from Edvardsson 
et al. (1993). Note, that none of the stars are significantly 
affected by interstellar reddening according to the color 
excesses derived from the Hp index and b ~ y. Accord- 
ing to the kinematical parameters of the stars as given in 
Edvardsson et al. (1993) they belong to either the thick 
disk or the old thin disk. The a-elements, e.g. O and Mg, 
are somewhat enhanced in the stars, ranging from [a/Fe] 
~ 0.15 in HR8181, HD 68284 and HD 130551 to [a/Fe] 
~ 0.25 in HR2883 and HR3578. 



3. Observations and data reduction 

The observations were carried out with the ESO Coude 
Echelle Spectrometer (CES) in three different periods: 
October 22-27, 1992, June 6-8, 1993, and February 5- 
9, 1995. In Oct. 92 and Feb. 95 the CAT 1.4m telescope 
was applied, whereas in June 93 the 3.6m telescope was 
feeding the CES throu gh a 3 5m long fiber and an image 
slicer (D'Odorico et al. 1989). On all occasions the detec- 
tor was a front illuminated Ford Aerospace 2048 x 2048 
CCD with 15 /im pixels. The CES camera has a dispersion 



of 1.88 Amm~^ at 6700 A providing a spectral coverage 
of about 58 A and 0.0282 Apixel"^. Different settings for 
the central wavelength were applied on the various nights 
in order to avoid any signatures in the spectrum resulting 
from a possible improper flatfielding. 

The main targets for the June 93 observations were 
some halo turnoff stars, for which the results have been 



presented by Smith et al. (1998) including a description 
of the reductions of the image slicer observations. The 
Oct. 92 and Feb. 95 observations have been reduced in 
much the same way. Subtraction of background, flat field 
correction and extraction of spectra were performed with 
standard tasks in IRAF. The wavelength calibration was 
based on 28 thorium lines well distributed over the wave- 
length region (6676 - 6734 A). A second-order polyno- 
mial was adopted for the dispersion solution resulting in 
a typical rms deviation from the fit of 0.0015 A. Further- 
more, the FWHM of the thorium lines were measured and 
found to vary by less than 4% along the spectrum. This 
near-constancy of the instrument profile is important in 
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Fig. 1. The observed spectra in the 6702 ~ 6710 A region. 
Spectral windows used in defining the continuum and in 
the x^ analysis of the lines are marked with C and L, 
respectively. Line identifications are shown for HR 3578 



connection with the determination of instrumental and 
stellar line broadening from Fel lines near the Lii line. 
For the Oct. 92 and Feb. 95 observations the FWHM of 
the Th lines corresponds to a resolution of i? = 105 000, 



Table 2. Number of observations and exposure times for 
the program stars 



ID 


Oct. 92 


June 93 


Feb.95 




CAT 1.4m 


3.6m 


CAT 1.4m 


HR2883 


4 X 40 min 






HR3578 


3 X 60 min 


2 X 25 min 


2 X 45 min 


HR8181 


4 X 30 min 


2 X 10 min 




HD 68284 






8 X 60 min 


HD 130551 






6 X 50 min 



whereas the June 93 image sclicer spectra have a resolu- 
tion of i? = 115 000. Finally, the spectra were normalized 
to an approximate level of 1 by fitting a fifth-order cu- 
bic spline function to the continuum, and corrected for 
the radial velocity shift using the accurate wavelengths 
(6703.567 and 6705. 102A) given by Nave et al. ( |1995D for 
the two Fe i lines close to the the Li line. 

Table 2 lists the the number of observations and ex- 
posure times as distributed over the three periods. Note, 
that HR 3578 has been observed in all three periods and 
HR 8181 in two of them. The individual spectra show ex- 
cellent agreement and have therefore been co-added to ob- 
tain the final spectra shown in Fig. |]for the 6702 - 6710A 
region. These are the spectra used for the determination 
of the ^Li/^Li ratio, but in the case of HR 3578 we have 
also derived the ratio separately for each observing period 
in order to check for possible systematic differences. 

4. The isotope ratio ®Li/^Li 

4.I. The center- of- gravity wavelength of the Lil line 

As discussed by e.g. Smith et al. ( 199^ ) the ^Li/'^Li ra- 
tio can be determined by two methods: From the center- 
of-gravity {cog) of the Li i 6708A line or from a detailed 
model atmosphere synthesis of the profile. The isotopic 
shift of the ^Li doublet is -1-0.158 A relative to the ^Li dou- 
blet. Addition of ^Li therefore shifts the Lil line to longer 
wavelengths and increases the FWHM. The cog-method 
relies in principle on a very simple and straightforward 
measurement, but its accuracy is limited by possible er- 
rors in the laboratory wavelengths of the lithium line and 
the reference lines needed to correct for the radial velocity 
shift of the star. The errors in the laboratory wavelengths 
are typically ±2mA. Differences in convective blueshifts 



may well be of the same order of size (Dravins 1987). To 
this should be added the uncertainty in measuring the cog- 
X for an asymmetric lithium line, which is slightly blended 
by a weak Fel line in the blue wing. According to our 
experience the cog-X will be uncertain by about ±5 mA. 
This translates to a one-sigma error of ±0.03 in ^Li/^Li, 
which is inferior to what may be obtained with the pro- 
file method. Hence, only the profile method will be used, 
and the wavelength of the Li i line will be considered as a 
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Fig. 2. The solar flux spectrum around the Lii line. The 
full drawn line is the observed spectrum from the Solar 
Flux Atlas of Kuruez et al. (1984). The dashed line is the 



synthetic spectrum computed for the model atmosphere of 
the Sun. The dotted line is the synthetic spectrum without 
the Li I line 



free parameter in the comparison between synthetic and 
observed profiles. 

^.2. Synthetic spectra 

The synthetic spectra have been obtained with the Up- 
psala Synthetic Spectrum Package, which computes the 
flux in a given wavelength region for an OSMARC model 
atmosphere of the star (Edvardsson et al. 1993 ). LTE is 
assumed and thermal and microturbulent line broadening 
as well as pressure broadening is included. The resulting 
spectrum is folded with a line broadening function that 



can be either Gaussian, radial-tangential (Gray 1978), ro- 
tational (Gray 1992) or any combination of these func- 
tions. As a check, one star (HD 68284) has also been ana- 
lyzed with a Kuruez ATLAS9 model atmosphere and the 
SYNTHE code (Kuruez |l99|). 

Wavelengths and ^/-values for the ^Li and ^Li compo- 
nents of the Li I doublet are taken from Table 3 of Smith 
et al. ( |1998D . As discussed by Smith et al. these values 
are not a source of significant errors in connection with 
determinations of the isotopic abundances of Li. 

In the case of halo turnoff stars with [Fe/H] < —1.5 
there is no significant blending of the Li i line by other 
lines, but already at [Fe/H] ~ —1.0 one has to worry 
about contributions from other absorption lines. To get 
more detailed information about line blending we there- 
fore started by synthesizing the solar spectrum, which is 
well suited for studying this problem because of the weak- 
ness of the Li I line. Fig. shows the solar flux spectrum 
from the atlas of Kuruez et al. (1984) compared to a syn- 
thetic spectrum based on the OSMARC solar model and 
a list of CN and metal lines from Miiller et al. (1975), who 
obtained loge(Li) = 1.0 and ^Li/^Li ~ 0.0 from a synthe- 
sis analysis of the solar Li i feature. We adopted these val- 
ues and adjusted the gf-values of the other lines to get the 



best possible fit of the synthetic spectrum (broadened by 
a radial-tangential function with a FWHM = 2.5kms~^) 
to the Solar Flux Atlas. As seen from Fig. the fit is 
quite satisfactory. The main problem is the two uniden- 
tified weak lines at 6708.02 and 6708.28 A, which have 
equivalent widths of 0.6 and 1.1 mA, respectively. The first 
line has nearly the same wavelength as the weak compo- 
nent of the ^Li doublet and therefore makes the determi- 
nation of ®Li/^Li very tricky for solar-type metallicities, 
whereas the unidentified line at 6708.28 A is further away 
and only affects the determination of the lithium isotope 
ratio marginally. 

Assuming that the two unidentified lines are neutral 
metal lines they will have equivalent widths less than 0.15 
and 0.3 mA, respectively, in solar-type stars with [Fe/H] ^ 
—0.6. This is too small to have any significant effect on 
the determination of the lithium isotope ratio. If the lines 
are ionized metal lines they will be somewhat stronger 
and affect the determination of ^Li/^Li marginally at the 
1% level. The probability for lines in solar-type spectra 
being from ionized metals instead of neutral is, however, 
small, and we have therefore chosen to exclude the two 
unidentified lines from our model atmosphere synthesis of 
the metal-poor disk stars. 

The CN lines seen in Fig. play no role in the spec- 
tra of the program stars mainly because both C and N 
scales almost linearly with Fe and partly because the 
stars have somewhat higher effective temperatures than 
the Sun. Hence, there is no reason to include these lines. 
The only line that has a significant effect on the synthesis 
of the lithium line is the Fei line at 6707.43 A. This line 
was included with the g/-value derived from the fit to the 
solar spectrum. 

As discussed in Sect. 4.3 the Fei lines at 6703.6 and 
6705.1 A have been used to determine the instrumental 
and stellar atmospheric line broadening. The two lines are 
blended by several faint CN lines in the solar spectrum 



as seen from Fig. 2 of Brault & Miiller (|1975|), but the 
synthesis of the solar and stellar spectra shows that these 
lines disappear beyond detection in the program stars. The 
same technique has been used to define a number of spec- 
tral windows practically free of lines in the program stars. 
These regions are marked by 'C in Fig. ^ and are used to 
set the continuum and to estimate the S/N in connection 
with the chi-square analysis of the lines. 

4-3. x^ analysis 

The synthetic spectra were computed for model atmo- 
spheres with the parameters given in Table 1. Line broad- 
ening due to macroturbulence and rotation was deter- 
mined from the two Fei lines at 6703.6 and 6705.1 A with 
Xexc — 2.76 and 4.61 eV, respectively. The synthetic spec- 
trum was first folded with a Gaussian function represent- 
ing the instrumental profile and then with either a Gaus- 
sian, a radial-tangential or a rotational profile. Various 
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5.8 5.2 6.6 7.0 

Tg (km/s) 

Fig. 3. Variation of the x^ of the fit to the Fei 6705.1 A 
line as a function of the FWHM of the Gaussian broaden- 
ing function 

combinations of these profiles was also tried, but in no 
cases the fit to the iron lines was better than that obtained 
with a Gaussian function. The final analysis was there- 
fore carried out with a single Gaussian representing the 
combined effect of instrumental and stellar atmospheric 
broadening. 

The FWHM of the Gaussian broadening profile, Tg, 
has been determined with the following procedure. First 
the continuum is set from the two windows on each side 
of a Fe I line and at the same time the standard deviation, 
a — {S/N)^^ of the spectrum is estimated. Then the 
chi-square function is computed: 

^(0^-^i)^ 



X 



n- 



(1) 



where 0^ is the observed spectral point and Si is the 
synthesis. The summation is performed over the spectral 
interval marked by 'L' in Fig. 1 corresponding to N = 25 
datapoints. In addition to Tq, the equivalent width and 
the exact wavelength of the Fe I line are considered as free 
parameters. Vq is varied in steps of 0.1 or 0.2 kms~^ and 
the other two parameters are optimized for each value of 
Tg to find the lowest x^- This results in a parabolic vari- 
ation of x^ BiS shown in Fig. ||. The most probable value 
of Tg corresponds to the minimum of x^i and Ax^ = 1,4 
and 9 correspond to the 1-, 2-, and 3-a confidence limits 



of determining Tg alone (Bevington & Robinson 1992). 

Table 3 summarizes the results of the x^ fitting of the 
two Fel lines. Note, that the reduced chi-square (Xred = 
X^ jv., where j/ = 22 is the number of degrees of freedom in 
the fit) is satisfactorily close to 1. Furthermore, the values 
of Tg from the two Fe i lines agree rather well, although 
there is a tendency that the 6705.1 A line gives slightly 
higher values of Tg than the 6703.6 A fine. 

Adopting a weighted average of Tg from Table 3 the 
lithium isotope ratio is determined from a x^ analysis 
of the fit between the computed and observed profile of 



Table 3. Results from the x^ analysis of the Fel lines 
at 6703.6 and 6705.1 A. Tg (6703) and Tg (6705) are the 
FWHM of the Gaussian broadening function (instrumen- 
tal -|- rotation -I- macroturbulence) applied to the syn- 
thetic lines. The errors given are the formal onc-sigma 
errors resulting from the x^ analysis 



ID 



S/N 



Tg (6703) 

[kms^^] 



2 
Xred 



Vg (6705) 
[kms^^] 



2 
Ared 



HR2883 


700 


5.70 ±.17 


0.78 


5.60 ±.11 


1.36 


HR3578 


1250 


5.35 ±.09 


1.07 


5.40 ±.06 


1.05 


HR8181 


1370 


5.25 ±.08 


1.17 


5.65 ±.05 


1.22 


HD 68284 


400 


5.40 ±.20 


0.96 


5.80 ±.14 


0.93 



HD 130551 830 6.30 ±.15 1.05 6.55 ±.10 0.96 




-0.04 -0.02 0.00 0.02 0.04 0.06 0.08 0.10 

f(6Li) 

Fig. 4. Variation of the ^ of the fit to the Lii 6707.8 A 
line as a function of the relative abundance of ®Li 



the Li I line. The free parameters in the fit are the total 
lithium abundance of the star, loge(Li), and the ^Li frac- 
tion, /(^Li) = A^(^Li)/iV(Li). Furthermore, a wavelength 
shift, AA, of the Li line relative to the wavelengths given 
in Table 3 of Smith et al. ( |1998| ) is allowed. The contin- 
uum and the S/N are determined from the two adjacent 
windows shown in Fig. n^ (resulting in nearly the same S/N 
values as given in Table 3) and the x^ analysis is carried 
out over the region marked by 'L', i.e. over N=33 data- 
points. Note, that the weak Fei line at 6707.43 A has only 
a small effect on the flux in this line region. 

The relative ^Li abundance is varied in steps of 0.01 
and the other two parameters, loge(Li) and AA, are op- 
timized for each value of /(^Li) to find the lowest x^- In 
order to study the behaviour of the x^ function around 
/(®Li) = 0.00 'negative' values of /(^Li) have been simu- 
lated by subtracting the line absorption coefficient due to 
^Li from the continuous absorption coefficient instead of 
adding it. 

The results of the x^ analysis are summarized in Table 
4 and the variation of x^ with /(^Li) is shown in Fig. ^ 
As seen, the three 'HR' stars have /(^Li) close to zero. 
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Fig. 5. The model atmosphere synthesis of the Fei 
6705. lA and the Lii 6707. SA line in the spectrum of 
HR8181. The datapoints are shown with open circles. In 
the upper figure the full drawn line corresponds to a Gaus- 
sian broadening parameter of Tg — 5.6 km s^^, whereas 
the dotted and dashed lines correspond to Tq = 5.0 
and 6.2knis^^, respectively. In the lower figure Tq = 
5.5 kms^^ has been applied. The full drawn line corre- 
sponds to /(®Li) = 0.0 and the dashed line to /(^Li) = 
0.05. Note, that when Tq and /(^Li) are varied the other 
free parameters in the fits, the wavelengths and the equiv- 
alent widths of the lines, have been optimized to get the 
best possible fits 



whereas ^Li has been detected in the two 'HD' stars at 
a high confidence level. As a further illustration of this 
result Figs. || and ^ show the fits to the Fe i and the Li i for 
two stars - one with /(^Li) ~ 0.0 and the other one with 
/(^Li) ~ 0.06 - and Fig. M shows a plot of the residuals 
in the observations after subtracting the ^Li and Fe i part 
of the synthesis. Although there is a disturbing periodic 
noise in the residuals with an amplitude of 0.1 to 0.2% in 
addition to the shot noise (the reason for which remains 
unexplained) a clear residual absorption at the wavelength 
of the ^Li doublet is seen in the spectra of HD 68284 and 
HD 130551. Note, that ahhough the S/N of the spectrum 



Fig. 6. Same as Fig. H for HD 130551. In the upper figure 
the full drawn line corresponds to a Gaussian broadening 
parameter of Fc — 6.6 km s^^, and the dotted and dashed 
lines to Fg — 6.0 and 7.2 kms^^, respectively. In the lower 
figure Tg = 6.6knis~^ has been applied. Here the full 
drawn line corresponds to /(^Li) = 0.06 and the dotted 
and dashed fines to /(^Li) = 0.00 and 0.10, respectively 



of HD 68284 is inferior to that of HD 130551, the error of 
/(^Li) is nearly the same for the two stars, because the Li 
line is about a factor of two stronger in the spectrum of 
HD 68284 than in the case of HD 130551. 

In order to check that there are no significant sys- 
tematic differences between the results obtained for the 
three observing periods the individual spectra of HR3578 
were analyzed separately including the determination of 
Tg. The x^ analysis gave the following values: /(^Li) = 
-0.018 ± 0.012, +0.008 ± 0.010, and -0.006 ± 0.014 for 
the Oct. 92, June 93 and Feb. 95 spectra. Within the 
errors these values agree satisfactorily with the value 
(/(^Li) = 0.000 ± 0.006) derived for the averaged spec- 
trum. 
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Table 4. Results from the y^ analysis of the Lii line. 
Yq is the weighted average of the FWHM of the Gaus- 
sian broadening function determined from the Fe i lines 
at 6703.6 and 6705.1 A (Table 3). loge(Li) is the total Li 
abundance and /(^Li) is the relative abundance of ^Li. 
The errors given are the formal one-sigma errors resulting 
from the x^ analysis. The last column gives the reduced 
chi-square of the fit 



ID Vg loge(Li) /C'Li) 

[kms-i] 



Xred 



HR2883 


5.7 


2.34 


0.010 ±.007 


0.99 


HR3578 


5.4 


2.05 


0.000 ±.006 


1.03 


HR8181 


5.5 


2.38 


0.004 ±.004 


0.89 


HD 68284 


5.6 


2.35 


0.041 ±.010 


1.09 


HD 130551 


6.5 


2.30 


0.056 ±.011 


1.22 



J^.J^. Systematic errors o/^Li/^Li 

The errors of /(^Li) given in Table 4 are purely statistical. 
In addition we must consider possible systematic errors 
resulting from approximations in the model atmospheres 
and the spectrum synthesis of the Fei and Lii lines. 

First we note that the exact values of the atmospheric 
parameters of the models are not critical for the lithium 
isotope ratio derived. An increase of Toff with +100 K 
changes loge(Li) with +0.07 dex but /(^Li) is practically 
unchanged. Changes in the gravity within reasonable lim- 
its have completely negligible effects. A decrease of the 
microturbulence by say 0.4kms~^ is compensated by a 
slight increase of Tq and the change of /(^Li) is less than 
0.003. A more significant, but still rather small change, 
results from the use of Kurucz's model atmospheres (Ku- 
rucz 1993) instead of the OSMARC models. As an exam- 



ple a model atmosphere with the parameters of HD 68284 
was interpolated between the closest models in the Ku- 
rucz ATLAS9 grid (re-computed with the "approximate" 
overshooting option switched off) and the SYNTHE code 
was used to compute profiles of the Fe i and Li i lines. This 
resulted in an increase of Fq from 5.6 to 5.9kms~^ and a 
change of /(^Li) from 0.041 to 0.033. Hence, the class of 
plane-parallel model atmospheres adopted has some effect 
on the value of the lithium isotope ratio derived, but not 
large enough to question the detection of ^Li in HD 68284 
and HD 130551. 

A more critical problem is whether classical plane- 
parallel model atmospheres are good enough for a deter- 
mination of /(^Li) with an accuracy of say ±0.01. It is 
well known that convective motions in real atmospheres 
produce slightly asymmetric line profiles, i.e. curved bi- 
sectors with a shape that depends on the depth of line 
formation (e.g. Dravins 1987). The question is therefore 



if the atmospheric velocity broadening determined from 
the two Fei lines is also valid for the Lii line. Due to the 
low excitation potential of the lithium resonance line it 
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Fig. 7. The residuals of the observations after subtraction 
of the ^Li and Fei 6707.43 A part of the synthesis of the 
Lii line. For comparison the synthesis of the ^Li doublet 
is shown with a full drawn line 



is probably formed somewhat higher in the atmosphere 
than the iron lines. In the case of HD 130551 a drastic 
increase of Tq to 7.5kms~^ instead of the value derived 
from the Fei lines, (6.5kms~^) decreases /(^Li) to about 
zero. One could also imagine that the Li i line in HD 68284 
and HD 130551 has a red asymmetry that mimics the ^Li 
doublet although no such asymmetry is seen in the pro- 
files of the Fe i lines. However, such possible effects have to 
occur for HD 68284 and HD 130551 only, because we can 
not allow a reduction of /(^Li) for the other three stars 
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Table 5. Radial velocities of the stars as derived from 
the accurate wavelengths of the Fei hnes (6703.567 and 
6705. IO2A) measured by Nave et al. ( [1995D . The wave- 
lengths shifts given are obtained from the y^ fits of the 
synthetic spectra to the individual lines 



ID 



RV 



[kr 



AA(6703) AA(6705) AA(Lii) 



lAl 



.Al 



rA] 



HR2883 


-1-54.8 


-0.1 


+0.0 


4.3 


HR3578 


+119.4 


+1.3 


-1.4 


9.1 


HR8181 


-28.8 


+2.6 


-2.5 


6.7 


HD 68284 


+62.8 


+1.3 


-1.4 


1.8 


HD 130551 


+33.7 


+1.4 


-1.4 


0.1 



which have /(^Li) ~ 0.00, when the atmospheric velocity 
broadening of the Fe i lines is adopted. 

As seen from Table 1 there are indeed significant dif- 
ferences in the atmospheric parameters of the two 'HD' 
stars and the three 'HR' stars that could induce some dif- 
ferential effects in the broadening of the Fe i and Li i lines. 
HD 68284 has a lower gravity and HD 130551 has a higher 
Teff than the other three stars. A hint that there may be 
systematic differences in the convective pattern between 
the two groups of stars comes from the small changes of 
the laboratory wavelengths needed to optimize the x^ fit 
of the Fei and Lii lines (see Table 5). The apparent he- 
liocentric radial velocities (including gravitational redshift 
and convective blueshift) is determined from the two Fe i 
lines. Hence, the sum of the wavelength shifts for these two 
lines is zero by definition. As seen, the wavelength shift of 
the Li I line is always positive. The average value is 4.4 niA 
(corresponding to a redshift of 0.20 kms^^ of the Lii line 
relative to the Fei lines) and the rms scatter is 3.4mA. 
This is more than the expected error of the laboratory 
wavelengths. In particular, we note that the redshift for 
HD 68284 and HD 130551 are lower than the redshift for 
the other three stars. Although this could be accidental, 
it is a warning that there may be differential effects in the 
convective line broadening. Clearly, this problem should 
be further studied by applying recently constructed inho- 
mogeneous 3D hydrodynamical model atmospheres (As- 
plund et al. 1999) in the analysis of the Lii line. 



5. Discussion 

5.1. Stellar masses 

According to standard stellar models the depletion of 
lithium is a strong function of stellar mass (Pinsonneault 
et al. 1992). Hence, it is of considerable interest to deter- 



mine the masses of the stars. This can be done by com- 
paring Teff and absolute magnitude, My, with mass tracks 
from stellar evolution calculations. 

Using the apparent magnitudes given in Table 1 and 
parallaxes from The Hipparcos and Tycho Catalogues 
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Fig. 8. Position of the stars in the logToff-My diagram 
compared with mass tracks from VandenBerg et al. (1999). 
Masses are given in units of the solar mass. The full 
drawn Hnes refer to [Fe/H] = —0.61 (equal to the metal- 
licity of HD 68284 and HD 130551) and the dotted hnes 
to [Fe/H] — —0.71 (approximately equal to [Fe/H] of the 
other stars). Both sets of tracks have [a/Fe] — 0.3. The er- 
ror bars in the x-direction correspond to a{Tcff) — ±70 K 
and those in the y-direction to the errors of the Hipparcos 
parallaxes 

Table 6. Absolute magnitudes computed from my and 
Hipparcos parallaxes and masses and ages derived from 
stellar evolu tionary tracks and isochrones of VandenBerg 
et al. (|1999| ) 



ID 



Mv 



M/Mq Age [Gyr] 



HR2883 3.55 +.06 

HR3578 4.16 +.04 

HR8181 4.41 ±.01 

HD 68284 3.41 +.19 



1.02 +.02 7.6 ±0.6 

0.88 +.02 11.7 ±1.2 

0.96 ±.02 6.3 ±2.0 

1.07 ±.04 6.8 ±1.4 



HD 130551 3.77 ±.09 1.06 ±.02 



6.0 ±0.9 



(ESA 1997[ ) the absolute magnitudes are calculated (Ta- 
ble 6), and the stars are plotted in the logTeS-My dia- 
gram (Fig. g). The mass tracks shown are from the new, 
g-ele ment e nhanced, evolutionary models of VandenBerg 
et al. ( 1999|) . Interpolation between the mass tracks (tak- 
ing into account their dependence on [Fe/H]) leads to 
the masses given in Table 6, and from the corresponding 
isochrones the stellar ages given are obtained. 

The errors of the masses and ages given in Table 6 
are standard errors corresponding to the adopted errors 
of Toff and My- Additional errors may be present due 
to inadequate stellar models and uncertainties in the cal- 
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ibration of T^s and the bolometric correction. Such er- 
rors are, however, more systematic and are expected to 
affect all stars with about the same amount. Hence, we 
conclude from Table 6 that the two stars for which ^Li 
has been detected (HD 68284 and HD 130551) have sig- 
nificantly higher masses than the three stars with no ^Li 
present in their atmospheres. This makes sense, because 
the depth of the convection zone of a star on the main se- 
quence decreases rapidly as a function of increasing mass. 
Hence, according to standard stellar models without mix- 
ing, the depletion of ^Li is less severe in the more massive 
stars. In this connection we note that although HD 68284 
is already on the subgiant branch and the coolest of the 
stars, it has spent most of its life as a main sequence star 
at Tcff ~ 6300 K . 

5.2. Galactic evolution and stellar depletion o/^Li 

Interpretation of the novel result of this paper - the detec- 
tion and quantitative measurement of the ^Li abundance 
in two old metal-poor disk stars - is contingent on two 
factors: (i) the expected evolution of the interstellar ^Li 
abundance with metallicity, and (ii) the depletion of the 
stellar ^Li abundance by the convective mixing that oc- 
curs in the pre-main sequence phase, and the additional 
depletion occurring on the main sequence. 

As is all too well known, prediction of Li depletion by 
main sequence stars and subgiants is an imprecise art. 
Standard models by Pinsonneault et al. ( 1992 ) predict 
loss of lithium in the pre-main sequence phase and no 
subsequent loss for stars of the mass of our quintet. De- 
pletions for masses of up to 0.85A^q and metallicities 
corresponding to [Fe/H] = —2.6 and —1.6 are computed 
by Pinsonneault et al. For their 0.85A1q model, the pre- 
dicted ^Li-depletions are negligible and the ^Li-depletions 
are 0.3 dex at [Fe/H] = — 1.6 and by extrapolation less 
than 0.05 dex at [Fe/H] = -2.6. Extrapolation to [Fe/H] 
~ —0.7 is uncertain but these depletions decrease with 
increasing mass such that our stars might be anticipated 
to have lost little, if any, ^Li. Cayrel et al. (1999b) report 



calculations that essentially confirm the above pre-main 
sequence depletions but predict a substantial continuing 
depletion of ^Li on the main sequence. At 7V1 = 0.85A^q 
and [Fe/H] = —1.5, a total ^Li depletion of about 0.7 dex 
is predicted in contrast to the 0.3 dex expected by Pin- 
sonneault et al. Consideration of non-standard physics, 
especially rotationally-induced mixing will result in likely 
larger and as yet more uncertain depletions - see, for ex- 
ample, the state of the art calculations by Pinsonneault et 
al. In summary, depletion of ^Li is to be expected but, at 
present, the magnitude of this depletion is uncertain with 
even standard calculations unavailable for the mass and 
metallicity of our old disk stars. 

Encouraged by recent observations of ^Li, Be, and B 
several predictions about the galactic chemical evolution 
of Li, Be, and B have appeared. Behind such predictions 



are assumptions about the nucleosynthetic processes of Li, 
Be, and B manufacture that demand assumptions about 
the early Galaxy, especially about the cosmic rays that 
permeated the halo and then the disk. Qualitatively, the 
key nucleosynthetic processes are known: (i) the Big Bang 
provided only the ^Li (in addition to H, ^H, "^He, and 
''He) that is widely considered to account for the observed 
Li abundance of the warm halo stars, the so-called Spite 
plateau; (ii) interactions between standard GCR and ISM 
and/or interactions between fast C,N,0 nuclei from super- 
bubbles with H or He in the ISM provide Li, Be, and B by 
spallation processes (e.g., O + p ^ Be) and Li through the 
fusion process a + a — >^Li and ^Li; (iii) neutrino-induced 
spallation processes in Type II supernovae that may pro- 
vide ^Li and '^^B. 

A key facet of this suite of processes is that beryllium 
with ^Be as the single stable isotope is produced solely by 
spallation of C,N,0 in fiight or at rest. Hence, observed 
beryllium abundances may be used to calibrate the yields 
of cosmic ray spallation. This is especially useful now that 
there are extensive measurements of the Be abundance in 
disk and halo dwarf stars. The relative yields of light nu- 
clides, for example ^Li to ^Be, are essentially independent 
of the cosmic ray spectrum unless there is a large excess of 
low energy cosmic rays {E < 30 MeV nucleon"'^) with re- 
spect to higher energy particles. This happy circumstance 
arises because above the similar threshold energies for the 
different processes (e.g., ^Be fromp -I- O and ^"B also from 
p + O), the spallation cross-sections are almost energy in- 
dependent. That ratios of yields are independent of the 
form of the (high ) ene rgy spectrum was well illustrated 
by Ramaty et al. (1996). At energies around the threshold 
energies for the various processes, the relative yields are 
energy and composition dependent. Moreover, Li produc- 
tion occurs also through a + a fusion reactions that do 
not synthesize Be and B. 

The B/Be ratio of halo stars is consistent within mea- 
surement uncertainty with production by spallation: Dun- 
can et al. ( 19971) e stimated B/Be = 15 ± 3 and Garcia 
Lopez et al. ( |1998 ) from a similar dataset of HST spec- 
tra found B/Be = 17 ± 10. Relativistic cosmic rays and 
the suite of (p,a) on (C,N,0) proces ses ar e predicted to 
give B/Be ~ 14 - see Ramaty et al. ( 1996| ) for predicted 
B/Be ratios as a function of cosmic ray energy and com- 
position. It has long been known that spallation by rela- 
tivistic cosmic rays is inadequate to account for the solar 
system's '^^B/'^'^B ratio which at 4.05 exceeds the predic- 
tion of about 2.5. Low energy spallation or a contribution 
from Type II supernovae are needed to resolve this dis- 
crepancy. 

These uncertainties aside, the Be observations are a 
reasonably firm basis from which to predict the ^Li abun- 



dances provided by spallation. Smith et al. (1998) dis- 
cussed the prediction of ^Li abundances from observed 
Be abundances - see the long-dashed line in Fig. g, where 
beryllium abundances are taken from Gilmore et al. (1992) 
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and Boesgaard et al. ( |l999b[ ). Predicted ^Li abundances 
are about a factor of 10 less than the observed ^Li abun- 
dances in the halo stars HD 84937 and BD -h26°3578 but 
exceed the ^Li abundances reported here for the disk stars 
HD 68284 and HD 130551 by about a factor of 3. Since 
^Li has almost certainly been depleted during pre-main 
sequence evolution and possibly during residence on the 
main sequence, the initial or interstellar ^Li abundance for 
the halo stars was higher than now observed. The required 
additional ^Li is probably primarily a product of cosmic 
ray a + a fusion production. 

Predictions of the growth of ^Li in the Galaxy made 
recently by Fields & Olive (1999a,b), Vangioni-Flam et al. 
( |1999D and Yoshii et al. ( |1997| ) are shown in Fig. |. These 
call on the same production processes but in different pro- 
portions. 

Fields & Olive discuss what they term the standard 
picture of galactic cosmic ray nucleosynthesis in a model 
galaxy. The key assumptions are that the cosmic rays al- 
ways had the composition of the ambient interstellar gas 
(i.e., they were very CNO-poor early in the life of the 
Galaxy), the energy spectrum of the cosmic rays was that 
measured for contemporary cosmic rays in the solar neigh- 
borhood (i.e., relativistic energies are dominant), and the 
cosmic ray flux has been proportional to the local super- 
nova rate, and scaled so that solar abundances of ^Li, B, 
and ^''B are reproduced. These assumptions with a sim- 
ple chemical evolution code (Fields & Olive report results 
for the canonical closed box) lead to the predicted run 
of the ^Li abundances with [Fe/H] where iron is a prod- 
uct of stellar nucleosynthesis with yields from Woosley & 



Weaver (1995) and a standard initial mass function. The 



key novel ingredient in the otherwise familiar calculation 
is the incorporation of recent measurements of the oxygen 
abundance in halo stars that indicate [O/Fe] increasing 
with decreasing [Fe/H] (Israelian et al. 1998; Boesgaard 
et al. 1999a). A higher O abundance increases the yields 



of spallation products. (A 'fudge' is needed as the O/Fe 
from these recent observations is considerably higher than 
predicted for Type H supernovae). The predicted ^Li vs 
[Fe/H] relation is shown by the dotted line in Fig. |^. A 
large part of the increase in the ^Li prediction at low 
[Fe/H] over the simple expectation from spallation is due 
to the inclusion of the a -|- a fusion reactions but the use of 
the observed O abundances through the associated contri- 
bution from p + O spallation appears necessary to match 
the observed ^Li abundances of the halo stars. Fields & 
Olive adjust their model to reproduce the solar ^Li which 
also accounts well for the Be and B abundances of the sun, 
disk and halo stars. 

Vangioni-Flam et al. (1999) incorporate a different mix 
of the light element producing processes into their chemi- 
cal evolution model. In particular, they invoke low energy 
nuclei that they associate with the acceleration of super- 
novae ejecta in the superbubbles created collectively by 
winds from the massive stars in OB associations. ('Low 




[Fe/H] 

Fig. 9. Abundances of lithium and beryllium as a function 
of [Fe/H] for 9 halo stars from Smith et al. ( |1998| ) and 5 
disk stars from the present paper. Open circles indicate the 
total Li abundance and filled circles the ^Li abundance or 
an upper limit. Open squares refer to the Be abundance 



adopted from Gilmore et al. (1992) and Boesgaard et al. 
(19991). The big symbols indicate meteoritic abundances 
from Anders & Grevesse (1989). The upper full drawn 
line is a fit to the 'Spite plateau' of lithium abundances 
for [Fe/H] < —1.5 and to the upper envelope of the Li 
abundance distribution for disk stars (see fig. 7 of Lambert 
et al. 1991). The lower full drawn line is a linear fit to the 



beryllium abundances with a slope of one, and the long- 
dashed line shows the corresponding relation for ^Li if 
^Li/Be = 5.8 as found in meteorites and as predicted from 
spallation of CNO nuclei by high energy cosmic rays. The 
dotted line rep resents t he evolution of ^Li in the model of 
Fields & Olive ( 1999a, h| ), the d ashed-dotted line the model 
of Vangioni-Flam et al. ( |1999 ) and the short-dashed lines 
refer to the model of Yoshii et al. ( 1997 ) for the halo and 
the disk, respectively 



energy' refers to energies close to the threshold energies 
of the spallation and fusion reactions.) A key point about 
this component is that the He, C, and O abundances of 
the ejecta are considered to be much higher than in the 
halo interstellar medium and, then, the dominant spalla- 
tion process is between (say) O in the ejecta and protons in 
the interstellar gas whereas in the standard picture (Fields 
& Olive 1999a, b) the leading process is between protons 



in the cosmic rays and (say) O in the interstellar gas. In 
Fig. H, we show predictions (dashed-dotted line) from a 
model adjusted to fit the measured ^Li abundances of the 
two halo stars. This model predicts a ^Li abundance at 
[Fe/H] = that exceeds slightly the solar abundance. 
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The close correspondence between the two predictions 
is unlikely to be a fair measure of the uncertainties in 
predicting the ^Li abundance of 1 Mq, [Fe/H] ~ —0.6 
disk stars starting from either the solar ^Li abundance or 
the ^Li abundance of halo stars. While the ^Li contribu- 
tion from spallation by galactic high-energy cosmic rays 
is rather well constrained by the observed Be abundance, 
there are no comparable constraints on the contributions 
of the fusion reactions and of spallation by low energy 
cosmic rays. 

That the range of permissible predictions is wider than 
perhaps suggested by the above two recent papers issug- 
gested by an earlier discussion by Yoshii et al. (1997). 
The prediction shown in Fig. ^ is from their Fig. 2 [J for a 
model that considers high-energy cosmic rays with cosmic 
ray protons and alphas spallating interstellar C,N, and O 
nuclei as well as a -|- a reactions. The cosmic ray flux was 
assumed to increase with decreasing metallicity. Different 
models are adopted for the halo and disk. This model pre- 
dicts a rather shallow decline of the ^Li abundance in the 
halo, and a steeper increase in the disk. Almost all of the 
^Li in the halo is the product of the a + a reactions. The 
prediction fails by about 0.4 dex to account for the solar 
^Li abundance, so that the discrepancy between predic- 
tion and observation for our disk stars might be larger for 
a revised model that did reproduce the solar ^Li abun- 
dance. 

Our old disk stars with detectable ^Li have Li abun- 
dances slightly in excess of the Spite plateau. If, as stan- 
dard models of pre-main sequence and main sequence evo- 
lution predict, the depletion of ^Li has been extremely 
slight, we may use the observed abundance and the pre- 
diction that cosmic ray production of the Li isotopes gives 
an isotopic ratio ^Li/^Li ~ 1.5 to estimate the contribu- 
tion of ^Li from cosmic rays. Consider HD 68284 with a Li 
abundance loge(Li) = 2.35. This is higher than the Spite 
plateau of loge(Li) = 2.21 (Smith et al. 1998). On the 
assumption that plateau stars have not depleted Li, the 
increase of Li in HD 68284 corresponds to Li ~ 60 on the 
scale H =10^^. If cosmic rays were entirely responsible 
for this increase, a division ^Li ~ 36 and ^Li ~ 24 is ap- 
propriate for a production ratio ^Li/^Li ~ 1.5. The Be 
abundance implies ^Li ~ 30 so that at this metallicity 
spallation rather than fusion reactions may be dominant. 
The observed ^Li abundance is ^Li ~ 10. (HD 130551 pro- 
vides similar figures.) Given that ^Li has assuredly been 
depleted to at least a modest extent, this elementary dis- 
section of the observed Li abundance reveals no obvious 
difficulty with a cosmic ray contribution to the Li isotopes. 



Recently, Ryan et al. ( 199E ) have argued on the basis 
of lithium abundances of Spite plateau stars in the range 
—3.5 < [Fe/H] < —2.3 that the plateau has a metallicity 



^ In converting [O/H] to [Fe/H] we have adopted [O/H] = 
[Fe/H] + 0.5 for [Fe/H] < -1.0 and [O/H] = 0.5 [Fe/H] for 
[Fe/H] > -1.0 



dependence due to the manufacture of the Li isotopes by 
cosmic rays. They consider the primordial abundance to 
be loge(Li) ~ 2.00 at [Fe/H] = -3.5. If Smith et al.'s 
Toff-scale is adopted, this abundance is raised to be about 
2.12 according to 4 stars common to both analyses. This 
and Ryan et al.'s metallicity dependence predict HD 68284 
to have a Li abundance of 2.44 which is similar to the 
observed value of 2.35. Relative to a plateau of 2.12, the 
observed abundance implies cosmic rays have added Li in 
the proportion ^Li = 54 and ^Li = 36 which are consistent 
with our observations provided that ^Li has been depleted 
by about 0.6 dex. At some point in the evolution of the 
Galaxy, sources (presumably stellar) contributed ^Li with 
little or no ^Li in order to raise the ^Li/^Li ratio to the 
solar ratio of 12.5. Inclusion of such a contribution in the 
above argument reduces the ^Li inferred from the increase 
in Li abundance over the plateau's value. 

The preceding argument may be inverted: the pre- 
dicted growth of ^Li with [Fe/H] may be used to infer 
the ^Li abundance. For example, the models proposed 
by Fields & Olive, and Vangioni-Flam et al. predict a 
^Li abundance in the ISM at the birth of HD 68284 and 
HD 130551 of about 120, a factor of 12 greater than ob- 
served. Assuming again a production ratio of ^Li/^Li = 
1.5 the attendant ^Li abundance is 180. Added to the 
primordial ^Li abundance of 160 this implies a total ^Li 
abundance of 340 (loge(^Li) — 2.53), a value considerably 
greater than the observed value of about 200. The obvi- 
ous implications are that either the predicted ^Li abun- 
dance is greatly overestimated or ^Li has been depleted 
by about 0.2 dex. In sharp contrast, Yoshii et al.'s predic- 
tion is a ^Li abundance of about 30 providing a total ^Li 
abundance of 205 or loge(Li) = 2.31, in excellent agree- 
ment with the observed ^Li abundances of HD 68284 and 
HD 130551. This model fails, however, to account for the 
meteoritic ^Li abundance by a large amount. 

6. Concluding Remarks 

It is presently impossible to refine our conclusions because 
understanding of Li-depletion, especially of ®Li-depletion, 
is poor. This is unfortunate as there is a tantalizing hint 
from the few detections of ^Li with the support of the 
greater suite of upper limits to the ^Li abundance that 
the growth of the ^Li abundance in the Galaxy up to the 
metallicity of the old disk ([Fe/H] ~ —0.5) may have been 
very slight, a result in conflict with the most recent models 
of ^Li chemical evolution but in agreement with a model 
proposed by Yoshii et al. ( 1997 ). This tentative suggestion 
is dependent on the absence of severe ^Li-depletion in our 
disk stars. 

Further progress will need additional observations of 
^Li in halo and disk stars. Such observations for extremely 
metal-poor stars, say [Fe/H] < —3, will require a high- 
resolution spectrograph on an extremely large telescope. 
More common instrumentation will be adequate for ex- 



12 



P.E. Nissen et al.: Isotopic lithium abundances in five metal-poor disk stars 



panding the sample of more metal-rich halo and disk stars 
with a known ^Li abundance and, in this case, a very 
large sample of stars may be needed to map out the up- 
per envelope to the ^Li vs [Fe/H] relation which with the 
corresponding '''Li vs [Fe/H] envelope may suffice to de- 
termine empirically the evolution of the ''Li/^Li ratio. A 
by-product of the survey will be measures of the range of 
^Li-depletions experienced by low mass stars of differing 
metallicity that will serve as grist for the theoreticians' 
mills. 
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